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1.  Statement  of  the  Problem  Studied 

An  understanding  of  the  elementary  steps  in  the  decomposition  of  energetic  materials  is 
essential  to  prediction  of  their  performance,  long-term  stability,  and  sensitivity  to  heat 
and  shock.  A  variety  of  sophisticated  experimental  approaches  have  thus  been  employed 
over  the  years  for  these  purposes.  However,  even  for  the  best  available  experimental  and 
theoretical  methods,  gaining  this  fundamental  understanding  is  quite  challenging  owing 
to  the  size  and  complexity  of  the  molecules  and  the  extreme  conditions  under  which  these 
reactions  take  place.  Given  this  complexity,  at  every  step  we  must  be  guided  by  our 
chemical  intuition  in  choosing  the  experiments  we  will  pursue  and  in  interpreting  the 
results  we  obtain.  Recent  work  has  revealed  an  important  new  mechanism  in  molecular 
decomposition,  termed  “roaming,”  that  has  prompted  a  significant  shift  in  our  views  of 
the  dynamics  of  unimolecular  dissociation.1  9  This  development  shows  that  we  have 
overlooked  a  significant  pathway  in  chemical  reactivity,  and  it  compels  us  to  reframe 
how  we  see  the  decomposition  of  energized  molecules  in  general,  and  of  high-energy- 
density  materials  in  particular.  Our  goal  in  this  program  is  thus  to  explore  in  detail  the 
potentially  profound  impact  of  hitherto  overlooked  roaming  dynamics  in  decomposition 
of  energetic  materials.  This  is  achieved  using  a  suite  of  powerful,  complementary 
experimental  techniques. 

2.  Summary  of  the  Most  Important  Results 

2.1  Roaming-dynamics  and  electronic  relaxation  pathways  in  the  nitrate  radical 

The  nitrate  radical,  NO3,  is  an  extraordinary  molecule  that  has  fascinated  and 
challenged  chemists  for  many  years.10'11  It  possesses  several  low-lying  electronically 
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excited  states  showing  complex  vibronic  interactions  accessible  via  visible  excitation.  " 
Its  photochemistry  is  also  quite  unusual,  with  two  dissociation  channels:  NO  +  O2  and 
NO2  +  O,  showing  nearly  the  same  energy  threshold  and  accessible  in  the  visible16,17'18 , 
giving  it  an  important  role  in  the  atmosphere.19  In  hindsight,  the  close  connection 
between  these  two  reaction  channels  may  be  seen  as  a  portent  of  the  recent  recognition  of 
the  key  role  of  roaming  dynamics  in  this  system.  ’  Roaming  in  the  case  of  NO3  is  quite 
unusual,  however.  In  some  ways  it  is  analogous  to  nitromethane,  in  that  it  involves  a 
roaming  isomerization  process.  However,  NO3  clearly  has  been  shown  to  involve  excited 
electronic  states  as  illustrated  in  Fig.  1.  In  recent  imaging  studies,  Grubb  et  al.  found 
that  rotationally- selected  NO  translational  energy  distributions  fell  into  two  distinct 
groups  associated  with  distinct  O2  vibrational  distributions.  The  low  rotational  levels 


were  formed  with  vibrationally  excited  O2,  while  higher  rotational  levels  were  associated 
with  vibrationally  colder  O?.  Based  on  a  careful  search  of  the  stationary  points  and  the 
intersections  among  the  relevant  potential  energy  surfaces,  Morokuma  and  co workers22 
suggested  the  following  picture  of  its  dissociation  dynamics:  After  optical  excitation  to 
the  D3  state  (the  lower  member  of  the  Jahn-Teller  distorted  E'  pair),  the  system  relaxes 
via  the  D3/D2  and  D2/D1  conical  intersections  (CIs)  to  arrive  on  the  dark  Di  state.  They 
argued  that  the  vibrationally  colder  O2  product  arises  from  roaming  dynamics  on  Di  (via 
the  Di  roaming  saddle  point  (RSP)  given  in  Fig  1),  while  the  vibrationally  excited  O2 
arises  via  dissociation  from  the  ground  electronic  state  via  the  Do  RSP.  Strong  support  for 
this  assignment  was  seen  in  the  lambda  doublet  propensities,  which  were  found  to  be 
opposite  and  perhaps  near-limiting  for  the  two  distributions,-  ’  _  as  well  as  in  subsequent 
quasiclassical  trajectory  (QCT)  calculations  on  separate,  uncoupled  D0  and  Di  potential 
surfaces.'  These  remarkable  results  prompt  several  questions:  What  is  the  branching  for 
dissociation  on  these  two  surfaces?  What  determines  this  branching?  Does  roaming  take 
place  exclusively  on  the  dark  state,  or  on  the  ground  state,  with  branching  upon  exit,  or 
on  both?  Is  the  branching  determined  at  the  intersection  between  the  Dj  and  Do  surfaces 
in  the  course  of  electronic  relaxation,  or  possibly  from  Do  up  to  D|  after  an  initial 
relaxation  directly  from  the  bright  D3  state  to  Do?  To  begin  addressing  these  questions, 
we  combined  state- selected  DC  slice  imaging  of  the  NO  product  following  visible  or  a 
combination  of  visible  and  infrared  excitation  from  a  CO2  laser,  along  with  additional 
QCT  calculations. 

Ideally,  to  probe  the  question  whether  roaming  takes  place  exclusively  on  the 
excited  state  or  the  ground  state,  one  might  contrast  dissociation  following  electronic 
excitation  with  that  arising  from  infrared  multiphoton  excitation.  Unfortunately,  NO3  has 
no  infrared  active  bands  that  overlap  the  CO2  laser  lines,  so  direct  CO2  laser  excitation  is 
not  readily  achieved.  However,  several  of  the  unusual  properties  of  NO3  allow  us  to  begin 
an  investigation  along  these  lines  using  state- selected  DC  slice  imaging  as  in  the  earlier 
experiments,  but  contrasting  pure  electronic  excitation  with  a  combination  of  electronic 
and  infrared  excitation.  We  obtained  the  total  translational  energy  distributions  (P(E)’s) 
from  the  images  and  although  we  see  distinct  differences  for  each  probed  rotational  level 
and  detection  pathway,  the  results  for  prompt  dissociation  with  purely  electronic 
excitation  show  little  difference  when  compared  to  dissociation  with  a  combination  of 
visible  and  infrared  after  a  0.7  (is  delay.  We  first  obtain  an  estimate  of  the  branching 
between  the  two  dissociation  paths  based  on  the  measured  distributions.  Our  strategy  here 
is  first  to  create  composite  translational  energy  distributions  by  weighting  the 
distributions  using  the  appropriate  linestrength  factors  ,  then  fitting  the  composite 
distribution  with  separate  Do  and  Di  components  obtained  from  the  QCT  calculations 
from  Bowman  and  coworkers.  The  agreement  is  generally  satisfactory.  The  relative 
yields  for  Do:Di  from  the  fits  were  found  to  be  0.92,  0.49,  0.31  and  0.44  for  rotational 
levels  21.5,  30.5,  33.5  and  36.5,  respectively.  Lower  rotational  levels  are  associated 
exclusively  with  Do.  This  state-dependent  branching  was  then  used,  with  linear 
interpolation  over  the  rotational  distributions  given  by  Wittig  and  coworkers17,  to  obtain  a 
final  branching  for  dissociation  of  6.1:1  Do:Di. 

We  now  turn  to  a  consideration  of  the  electronic  relaxation  and  dissociation 
dynamics.  The  main  possibilities,  highlighted  in  Fig.  1,  include:  initial  excitation  to  the 
D3  bright  state,  followed  by  internal  conversion  (IC1)  through  D2  and  Di  to  Do  via  the 
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conical  intersections  that  were  found  by  Morokuma  and  coworkers.'  This  pathway  is 
shown  as  “IC1”  in  Fig.  1.  The  branching  for  dissociation  along  D|  or  Do  that  we  measure 
may  then  reflect  the  branching  at  the  D|/D0  Cl  as  it  is  encountered  during  the  relaxation. 
This  is  the  pathway  implicit  in  the  treatment  of  Grubb  et  al.  Alternatively,  excitation  to 
the  bright  state  maybe  followed  by  IC  to  the  ground  state  directly  from  D3,  shown  as 
“ICO”  on  Fig.  1.  The  distinction  here  is  that  we  recognize  that  D3  and  Do  are  directly 

vibronically  coupled13,  so  the 
ground  state  need  not  be  accessed 
primarily  via  Di.  In  this  scenario, 
after  IC  to  Do,  the  system  may 
proceed  to  the  roaming  region, 
then  “roam  upstairs”  via  the  D0/Di 
Cl,  which  is  located  in  the  roaming 
region.  From  there  it  may  roam  on 
Dj  or  simply  dissociate.  If  the 
energy  is  above  the  O  atom  loss 
threshold,  decay  via  that  path  is 
strongly  preferred  both  for  Dj  and 
Do,  based  on  both  experiment16  and 
on  the  trajectory  calculations.24 

Taking  all  these  points  into 
consideration,  we  believe  the 
evidence  supports  the  second 
alternative  relaxation  pathway  mentioned  above,  i.e.,  direct  IC  from  the  bright  state  to  Do, 
with  Di  subsequently  accessed  from  the  ground  state,  probably  via  mixing  in  the  roaming 
region  where  these  surfaces  are  strongly  coupled  as  shown  in  Fig.  5.  The  points  in 
support  of  this  view  are  as  follows:  The  D3  and  Do  surfaces  possess  strong  vibronic 
coupling,  directly  impacting  both  the  ground  state  vibrational  level  structure  and  the 
bright  state  fluorescence  lifetimes,  while  that  between  the  bright  state  and  Di  is  very 
weak.14  Moreover,  the  fluorescence  lifetimes  suggest  a  dilution  of  the  bright  state  that  can 
only  be  accounted  for  by  the  ground  state  density  of  states.26'27  The  fact  that  our 
translational  energy  distributions  are  essentially  identical  whether  we  use  electronic 
excitation  with  probe  on  a  10  ns  timescale,  or  coupled  with  IR  excitation  0.7 
microseconds  later,  suggests  that  this  relaxation  is  complete  on  the  10  ns  timescale. 

2.3  IRMPD  imaging  of  nitroalkanes  and  alkyl  nitrites 

The  advent  of  infrared  multiphoton  dissociation  (IRMPD)  in  the  early  70’ s  has 
been  highly  instrumental  in  the  study  of  nascent  photofragments  originating  from 
dissociation  of  small  molecules  on  the  ground  electronic  state.  With  the  introduction  of 
powerful  TEA-CO2  lasers,  IRMPD  became  a  unique  tool  for  studying  unimolecular 
dissociation.  It  has  been  used  extensively  to  study  reactions  that  involve  concerted 
pathways  as  well  as  tight  transition  states  and  simple  bond  rupture.  IRMPD  has  more 
recently  been  exploited  for  mass  spectrometric  applications.  However,  roaming  dynamics 
are  likely  to  occur  from  the  ground  state  near  threshold,  so  IRMPD  offers  a  powerful  way 
to  probe  these  processes.  Roaming-mediated  isomerization  has  been  invoked  in 
dissociation  of  nitrobenzene15  and  nitromethane14'16  to  account  for  the  significant  yield  of 


Fig.  1.  Relevant  features  of  the  potential  energy 
surfaces  of  the  nitrate  radical. 
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the  NO  loss  channel  in  these  systems.  To-date,  however,  in  primary  photodissociation 
studies  of  roaming  dynamics,  only  UV  excitation  has  been  used  to  prepare  excited 
molecules  for  subsequent  investigation.  Since  reactions  involving  roaming  mechanisms 
largely  arise  from  the  ground  electronic  state,  IRMPD  studies  can  be  very  helpful  in 
providing  unambiguous  information  about  such  processes.  In  our  recent  studies,  we  have 
applied  IRMPD  with  powerful  state-resolved  ion  imaging  methods  for  the  first  time. 

We  have  recently  performed 
an  IRMPD/DC  slice  imaging  study  of 
nitroethane,  nitropropane, 

nitrobutane,  butylnitrite,  and 
isobutyl  nitrite  with  state-resolved 
probing  of  the  NO  product.  Infrared 
multiphoton  dissociation  (IRMPD), 
is  used  to  dissociate  the  molecules 
on  the  ground  electronic  state, 
which  is  a  significant  advantage  in 
that  we  can  be  confident  we  are 
studying  ground  state  processes, 
unlike  the  case  of  electronic 
excitation.  The  channels  that 
produced  NO  as  a  product  were 
studied  and  compared  with  each 
other  to  gain  insight  into  the 
isomerization  mechanism  and  the 
dissociation.  We  show  preliminary 
results  for  nitrobutane  and  butyl 
nitrite,  and  we  are  in  the  process  of 
analyzing  the  rich  datasets  we  have 
obtained  for  all  of  these  systems. 

Resonant  multiphoton 

ionization  (REMPI)  spectra  were  recorded  for  the  NO  product  of  each  system 
following  IRMPD  on  the  0-0  band  for  the  A(2£+)<-X(2nr)  transition  in  the  range  of 
44180  -  44450  cm4.  The  peaks  are  assigned  to  each  rotational  level  by  comparing 
to  a  simulated  NO  spectrum  obtained  using  LIFBASE.  We  show  the  results  for  butyl 
nitrite  and  nitrobutane  in  Fig.  2  by  way  of  example.  It  is  immediately  apparent  that 
the  butyl  nitrite  has  a  much  hotter  rotational  distribution  than  nitrobutane. 

The  total  translational  energy  distributions  were  also  obtained  from  the  DC 
slice  images  of  the  NO  product  from  the  IRMPD  of  these  nitroalkanes  and  alkyl 
nitrites.  The  images  were  recorded  for  the  nascent  NO  products  for  different  J  values 
from  the  dissociation.  Total  translational  energy  distributions  were  obtained  by 
using  the  conservation  of  mass  and  momentum.  In  figure  3  the  images  are  shown 
alongside  the  total  translational  energy  distribution  for  a  single  NO  quantum  state 
for  IRMPD  of  1-nitrobutane  and  butyl  nitrite.  The  obtained  total  translational 
energies  peak  at  zero,  which  suggest  barrierless  dissociation  of  the  molecules.  In  the 
images  we  see  a  significantly  greater  translational  energy  release  for  the  butyl 
nitrite  compared  to  nitrobutane,  consistent  with  the  REMPI  spectra.  Interestingly, 
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Fig.  2.  REMPI  spectra  for  NO  from  nitrobutane  and 
butyl  nitrite 
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these  results  are  exactly 
counter  to  the  trends  we 
observed  in  nitromethane 
and  methyl  nitrite. 


The 


average 


Fig.  3  DC  sliced  images  (right)  and  translational  energy 
distributions  (left)  for  NO  J=10.5  from  IRMPD  of  indicated 
molecule. 


translational  energy  release 
for  each  system  and  each 
rotational  level  was  obtained 
from  the  translational  energy 
distributions  and  these  are 
given  for  all  systems  studied 
and  all  rotational  levels,  along 
with  the  rotational 
temperatures,  in  Fig.  4.  The 
average  translational  energy 


release  for  nitroethane  ranges  from  4.9  kcal/mol  to  7.0  kcal/mol  (low-J  to  high-}) 
while  for  ethyl  nitrite  it  ranges  from  1.2  kcal/mol  to  1.4  kcal/mol.  The  limiting  value 
for  the  total  translational  energy  release  of  nitroethane  was  found  to  be  ~14 
kcal/mol  and  ~5  kcal/mol  for  ethyl  nitrite.  This  difference  in  the  translational 
energy  release  is  also  manifested  in  the  product  rotational  excitation:  1600K  for 
nitroethane  and  768  for  ethyl  nitrite.  We  note  also  that  for  nitroethane  (and  the 
other  nitro  systems),  there  is  a  correlation  between  average  translational  energy 
release  and  rotational  excitation  (Fig.  8),  consistent  with  an  impulsive  origin  of  the 
rotation.  Interestingly,  this  is  not  seen  for  the  nitrites.  These  results  show  the  same 
trends  as  the  IRMPD  results  of  nitromethane  and  methyl  nitrite  we  reported 
previously,  supporting  the  association  to  RMI [2].  The  isomerization  of 

nitroethane  to  ethyl  nitrite  has  to  compete  with  the  direct  C-N  bond  fission  channel 
and  an  alternative  pathway,  concerted  molecular  elimination  (CME)  to  HONO  + 
ethylene,  for  which  there  is  no  analogue  in  nitromethane.  Denis  et  al.  reported  a 
density  functional  study  of  the  stationary  points  relevant  to  2003.  They  found  the 
barrierless  simple  bond  fission  (SBF)  of  the  C-N  bond  at  53.8  kcal/mol  above  the 
nitroethane  minimum,  and  they  placed  the  TS  for  CME  about  12  kcal/mol  lower. 
They  found  a  tight  TS  for  isomerization  to  the  nitrite  6  kcal/mol  above  the  SBF 
channel  and  troubled  over  this,  as  they  also  found  a  similar  value  for  nitromethane 
but,  as  noted  by  Lin,  this  high  barrier  would  not  permit  effective  competition 
between  isomerization  and  SBF  as  had  been  seen  for  nitromethane  in  the  Wodtke 
and  Lee  experiments.  Denis  et  al.  concluded  in  any  case  that  the  CME  channel 
would  dominate  in  nitroethane  and  nitropropane,  and  took  solace  in  the  fact  that 
Wodtke  and  Lee  did  not  see  a  measurable  yield  of  NO  in  the  nitroethane  case. 
However,  our  results  clearly  show  that  isomerization  to  the  nitrite  does  occur,  as 
this  is  the  only  feasible  means  of  producing  NO  from  the  nitroethane.  Furthermore, 
the  likely  mechanism  is  RMI,  just  as  in  the  nitromethane  case.  However,  although 
the  branching  to  isomerization  is  substantial  for  nitromethane,  our  results  give  no 
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indication  as  to  the  branching  here:  it  could  be  and  likely  is  relatively  small.  Indeed, 
Wodtke  et  al.  reported  no  evidence  for  isomerization  in  IRMPD  of  nitroethane  or  2- 
nitropropane;  However  they  did  not  report  an  upper  bound  for  this  channel.  Given 
mass  spectrometric  interference  from  both  NO2  and  HONO  at  m/z=30  in  their  study, 
and  the  challenge  of  identifying  ions  associated  with  C2H5O,  it  is  reasonable  to 
imagine  that  branching  up  to  15%  could  easily  be  overlooked  in  those  experiments. 
In  any  case,  our  results  suggest  that  there  is  a  loose  roaming  transition  state 
analogous  to  that  found  by  Zhu  and  Lin  for  the  nitromethane  case,  or  else 
isomerization  would  be  very  unlikely,  as  suggested  by  Denis  et  al..  Nitroethane  and 
the  other  systems  studied  here  all  have  the  vicinal  carbon  that  allows  for  the  five- 
membered  transition  state  leading  to  CME.  In  a  sense  these  systems  have  aspects  in 
common  with  the  formaldehyde  potential  surface  in  which  roaming  was  first 
recognized.  In  formaldehyde,  the  SBF  threshold  and  the  roaming  TS  are  also  above  a 
CME  transition  state,  about  7  kcal/mol  higher  in  that  case.  The  extent  to  which 
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Fig.  4.  Average  translational  energy  distribution  change  with  rotational  level  and 
rotational  temperature. 


roaming  persists  above  the  SBF  threshold  is  strongly  system-dependent,  and  a 
subject  of  considerable  interest.  In  formaldehyde,  branching  to  the  roaming 
pathway  remains  significant  even  at  energies  10  kcal/mol  or  more  above  the 
roaming  threshold  [3  5].  In  the  nitrate  radical,  RMI  shuts  down  almost  immediately 
above  threshold.  In  nitromethane,  there  is  some  disagreement  among  the 
theoretical  treatments.  Homayoon  and  Bowman  used  full-dimensional  quasiclassical 
trajectory  calculations  and  saw  RMI  persist  more  than  5  kcal/mol  above  the  SBF 
threshold.  Annesley  et  al.  used  a  reduced  dimensionality  "rigid-body"  trajectory 
treatment  and  concluded  that  RMI  is  operative  only  up  to  2  kcal/mol  or  so  above 
threshold.They  also  reported  the  pressure  dependence  and  concluded  that  SBF  will 
dominate  at  high  pressures  (120  torr).  We  should  note  that,  following  isomerization 
to  the  nitrite,  pathways  leading  to  HNO  +  acetaldehyde  are  also  accessible,  either  via 
a  four-center  transition  state,  or  by  NO  roaming.  A  recent  report  of  ethyl  nitrite 
pyrolysis  using  chirped-pulse  mm-Wave  spectroscopy  and  kinetic  modeling  has 
concluded  that  roaming  dynamics  are  responsible  for  HNO  formation  in  that  case. 
Future  studies  to  determine  the  branching  accurately  for  these  systems  would  be  of 
considerable  interest. 


2.4  Chirped-pulse  probing  of  ethyl  nitrite  flash  pyrolysis 
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In  collaboration  with  the  Field  group  at  MIT,  we  have  begun  application  of 
the  powerful  technique  of  chirped-pulse  (CP)  spectroscopy  to  reaction  dynamics. 
We  illustrate  this  with  recent  results  obtained  at  MIT  using  CP  millimeter-wave 
spectroscopy  coupled  to  flash  pyrolysis  of  ethyl  nitrite.  This  approach  was  used  to 
measure  product  branching  and  extract  information  about  the  reaction  mechanism, 
guided  by  kinetic  modeling.  Pyrolysis  of  ethyl  nitrite,  CH3CH2ONO,  was  studied  in  a 
Chen  type  flash  pyrolysis  reactor  at  temperatures  of  1000-1800  K.  The  branching 
ratios  of  HNO,  CH2O,  and  CH3CHO  products  are  measured  and  compared  to  the 
kinetic  models  generated  by  the  Reaction  Mechanism  Generator  software.  We  find 
that  roaming  CH3CH2ONO  ->  CH3CHO  +  HNO  plays  an  important  role  in  the  thermal 
decomposition  of  ethyl  nitrite,  with  its  rate,  at  1000  K,  comparable  to  that  of  the 
radical  elimination  channel  CH3CH2ONO  ->  CH3CH2O  +  NO.  HNO  is  a  signature  of 
roaming  in  this  system. 


2.5  IRMPD  of  vinyl  chloride  yields  cold  vinylidene 

Velocity  map  imaging  of  the  infrared  multiphoton  dissociation  of  vinyl 
chloride  shows  the  formation  of  HC1  in  rotational  levels  below  J=10  that  are 
associated  with  the  3-center  elimination  pathway.  The  total  translational  energy 
release  was  observed  to  peak  at  3-5  kcal/mol,  which  is  consistent  with  the  low 
reverse  barrier  predicted  for  the  formation  of  HC1  with  vinylidene  co-products. 
Direct  dynamics  trajectory  studies  from  the  3-center  transition  state  reproduce  the 
observed  distributions  and  show  the  associated  vinylidene  is  formed  with  only 
modest  rotational  excitation,  precluding  Coriolis-induced  mixing  among  the  excited 
vibrational  levels  of  acetylene  that  would  lead  to  distribution  of  vinylidene 
character  into  many  vibrationally  mixed  acetylene  vibrational  levels.  The  results 
suggest  that  infrared  multiphoton  dissociation  of  vinyl  chloride  is  an  efficient  route 
to  synthesis  of  stable,  cold  vinylidene. 
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